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ABSTRACT 

We put constraints on the properties of the progenitors of peculiar calcium-rich tran- 
sients using the distribution of locations within their host galaxies. We confirm that 
this class of transients do not follow the galaxy stellar mass profile and are more likely 
to be found in remote locations of their apparent hosts. We test the hypothesis that 
these transients are from low metallicity progenitors by comparing their spatial dis- 
tributions with the predictions of self-consistent cosmological simulations that include 
star formation and chemical enrichment. We find that while metal-poor stars and 
our transient sample show a consistent preference for large offsets, metallicity alone 
cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps 
to improve the match, indicating these events may result from a very old metal-poor 
population. We also investigate the radial distribution of globular cluster systems, and 
show that they too are consistent with the class of calcium-rich transients. Because 
photometric upper limits exist for globular clusters for some members of the class, 
a production mechanism related to the dense environment of globular clusters is not 
favoured for the calcium-rich events. However the methods developed in this paper 
may be used in the future to constrain the effects of low metallicity on radially distant 
core-collapse events or help establish a correlation with globular clusters for other 
classes of peculiar explosions. 
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1 INTRODUCTION 

Supernovae (SNe) are luminous events caused by the explo- 
sions of stars. They are divided into observational subclasses 
according to line features from different chemical elements 
appearing in the spectra; type la supernovae (SNe la) are 
generally attributed to the thermonuclear disruptio n of a 
white dwarf (see e.g.. lHillebrandt fc Niemever 1120001 ). while 
type Ib/c supernovae (SNe Ib/c) and type II supernovae 
(SNe II) result from the gravitational collapse of the core 
of a mas sive star at th e end of its life (main sequence mass 
J>8 Mm . ISmarttl 120091 ) . However, details of the progenitor 
systems, the physics of the explosions, and the evolutionary 
paths leading to those explosions are still unclear and under 
investigation. 



E-mail: yuanfang@mso.anu.edu.au 



Automated wide field optical transient surv eys, e.g., 
the Texas Supernova Search (TSS, lQuimby||2006l) /RO TSE 
Supernova Verification Project (RSVP, I Yuan! I2010T). the 



Catalina Real-Time Transient Survey (CRTS, Drake et al 
2009), the Palomar Transient Factory CPTF. iRau et al 



2009) and the Panoramic Survey Tele scope and R apid Re- 
sponse System (Pan-STARRS, iKaiser et al.ll2002T ). have in 
recent years not only discovered hundreds of new SNe of 
all types, but discovered new types of transient events in 
hitherto unexplored regions of parameter space. 

An increasing number of unusual events are being dis- 
covered by these surveys, often occurring in host environ- 
ments poorly sampled by previous supernova surveys, which 
targeted nearby massive galaxies. For example, PTFlOops 
is a peculiar SN la ~ 148 kpc away from the nearest galaxy 
l|Maguire et al.ll2011I ). It is puzzling to also find core-collapse 
explosions at remote locations, given the expected low lo- 
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cal star-formation rate. The massive star progenitors of 
these core-colla pse events either have very low m etallicity 
(e.g., SN 1983K. lNiemela et al.|[l985l ; iPhillips et al.lll990l ) or 
have gained very h igh velocity and migrated outward from 
the galaxy centre (|Zinn et al.l I2011T ). Alternatively, a core- 
collapse could be delayed if it is the result of the merger of 
two white dwarfs (e.g., an ONeMg and a He white dwarf; cf. 
lYoon et al.l 120071 ). At present, it is not clear whether such 
mergers prefer old populations and/or low-metallicity envi- 
ronments; 

IPerets et al.l (|2010|) identified a group of faint events, 
spectroscopically similar to type Ib/c SNe and displaying 
prominent calcium lines in the nebular phase. Whether these 
events constitute a homogeneous class with a common phys- 
ical explosion mechanism remains unclear, as does the na- 
ture of that mechanism (thermonuclear vs. core-collapse). 
A core-collapse scenario has been put forward for one of 
these events, SN 2005cz, in which an 8 to 12 Mq pro- 
genitor lo st its hydrogen-r i ch en velope through binary in- 
teraction IjKawabata et al.l l2010f ). However, detailed stud- 
ies at and around the supernova locati on do not reveal ev - 
idence of a young stellar population (jPerets et al.l 120111 ). 
The core-collapse scenario for these transients also strug- 
gles to explain the large fraction found in early-type host 
galaxies. An alternative explanation for the small amount 
of 56 Ni synthesised in these explosions is a scenario involv- 
ing a h elium detonation on the surface of an acc r eting white 
dwarf llBildsten et al. 120071; IPerets et all lioiol ; IShen et all 
l2010l ; IWaldman et aUhoiltlSim et alj|2012l ). Both the core- 



collapse and the thermonuclear mod els can reproduce qual- 
itatively a large Ca/ O line ratio IjKawabata et alj 120101 : 
IWaldman et al.l 1201 ll ). but comprehensive spectral series 
modelling across the range of observed calcium-rich tran- 
sients is still lacking. 

Several more events were discovered by PTF that share 
some com mon properties with t his calcium-rich group. In 
particular, iKasliwal et al.l (|2012l ) identified a class of five 
objects (including two events from IPerets et al.ll2010l ) using 
both photometric and spectroscopic criteria. Interestingly, 
all these peculiar PTF transie nts are fou nd at large di stances 
from their host galaxies (Sullivan ct al. 2011; Maguir e et alj 

No single model so far has successfully reproduced 
all the observed features of the peculiar "calcium-rich" 
transients. In this paper, we take an alternative, model- 
independent approach to constrain the properties of their 
progenitors by statistically comparing their locations to the 
overall stellar population in their host galaxies. We define 
our sample selection in Section [2] We quantitatively eval- 
uate the transient locations relative to stellar mass distri- 
bution in Section [3] We compare the observed distribution 
with theoretical predictions of cosmological simulations in 
Section [4] and investigate the expected radial distribution of 
objects associated with globular clusters in Section [S] The 
implications are discussed in Section HJ] and we conclude in 
Section (7) 



IKasliwal et"afl l|2012l ). This will be referred to as the K12 
sample. This sample includ e s two events (SN 2005E and SN 
2007ke) from IPerets et alj (|2010l ). as well as three events 
selected from a sample of over a thousand transients discov- 
ered by PTF in its first two years of operation: PTF09dav, 
PTFlOiuv (aka SN 2010et), and PTFllbij. Compared to 
typical supernovae, the defining characteristics are a mod- 
erately faint peak luminosity (absolute R magnitude of 
~ —16), relatively fast rise/decay (~12 - 15 days), normal 
photospheric velocity, a quick spectral transition into the 
nebular phase, and nebular spectra dominated by calcium 
emission. 

Another five events (SN 2000ds, 2001co, 2003H, 2003dg 
and 2003dr) have been spectroscop ically selected as pecu- 
liar calcium - rich Type Ib/c SNe (|Filippenko et a.1.1 120031 : 
IPerets et a.1.1 |2010| ). Their photospheric spectra resemble 
that of typical SNe Ib/c, lacking silicon, hydrogen and 
sometimes helium lines. Their nebular spectra are domi- 
nated by calcium emission. While it is impossible to verify 
the photometric properties for these objects, there is evi- 
dence that they ar e sub-luminous compared to normal SNe 
|Perets et al.ll2010l ). and consistent with the objects in the 
K12 sample. We therefore include them in the "calcium- 
rich" group to maximise o ur sample size for sta tistical study. 

All seven events from IPerets et al.l (J2010I) were discov- 
ered/ob served by the Lick Ob servatory Supernova Search 
(LOSS, iFilippenko et al.ll200ll) . We regard them as a rea- 
sonably complete and uniform set, although we will discuss 
possible selection effects in both the PTF and LOSS sam- 
ples in Section lrTTl Observatio ns of SN 2005cz, a singl e event 
discovered by another group (|Kawabata et al.1120101 ). show 
it to be a convincing member of the "calcium-rich" class. It 
is, however, excluded in our study to avoid further compli- 
cating the und erstanding of its selec tion bias. 

We follow IKasliwal et alj (|2012l ) in identifying the host 
as the nearest galaxy in projected distance at a redshift con- 
sistent with the transient. The measured distance between 
a transient and the centre of its host is thus a lower limit on 
the real distance. 



3 LOCATIONS COMPARED TO STELLAR 

MASS 

In this section we compare the distribution of the "calcium- 
rich" transients with the stellar mass distribution of their 
host galaxies. We use A'-band photometry from the Two 
Micr on All Sky Survey (2MASS) catalog l|Skrutskie et all 
2006|) to trace the overall stellar mass. 

We analyse images from the 2MASS All-Sky Data Re- 
lease. For the host of SN 2005E and SN 2007ke, several over- 
lapping images must be stitched together to get good cov- 
erage of the target galaxies. We subtract bright stars near 
the galaxies using the IRAF0 DAOPHOT package. We ex- 
tract the moments of the galaxy isophotes using SExtractor 



2 SAMPLE OF TRANSIENTS 

We first consider the collection of five transients selected 
based on both photometric and spectroscopic properties by 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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IJBertin fc Arnoutslll996r ). then use the IRAF package EL- 
LIPSE to extract fluxes in elliptical annuli, with the semi- 
major axis increased in one-pixel increments and a common 
centroid, ellipticity, and position angle for all annuli fixed at 
the values obtained from SExtractor. We fit the smooth part 
of the surface brightness profile between one- half and one 
Kron radius (2.5 ti mes the first moment radius, |KronJ[l980|) 
with a Sersic model (Sersic 1968), but allow the Sersic index 
to vary freely. We integrate the resulting surface brightness 
profile to estimate the total flux of the galaxy. Fractions of 
stellar light enclosed by elliptical apertures defined by the 
transient locations are calculated and listed in Table [T] 

SN 2003H lies between a pair of interacting galaxies. 
We ignore the brightening on the outskirts of IC 2163 due to 
contamination from NGC 2207. Nevertheless, it is likely that 
the local environment has been affected by the interaction 
and therefore cannot be fully represented by its distance to 
either galaxy in the pair. 

For SN 2007ke, both the transient and the identified 
host MCG +07-07-003 lie in the bright galaxy NGC 1129. 
We remove the background light from NGC 1129 by mod- 
elling its surface brightness profile in the same way as de- 
scribed above. If NGC 1129 is the host galaxy, the location 
of SN 2007ke will not be as extreme (within ~50% of light), 
but uncertainty in this single event does not change our main 
conclusions. 

If the surface brightness in the host galaxy drops faster 
(slower) than our model in the outer regions of the galaxy, 
the fraction of light enclosed by the elliptical isophote pass- 
ing through the transient's location will be larger (smaller) 
than our calculation. We estimate the systematic uncer- 
tainty in the surface brightness by considering profiles with 
Sersic indices of 0.2 (steep profile) and 4 (flat profile) outside 
the fitted region. 

Figure[l]shows that the "calcium-rich" transients do not 
follow the distribution of A"-band light, a proxy for the stel- 
lar mass, in their host galaxies. Instead, they have a strong 
preference for remote locations. A Kolmogorov-Smirnov (K- 
S) test shows that even the limiting case of a flat profile 
(Sersic index of 4) is inconsistent at greater than 99% con- 
fidence with a scenario in which the transient locations uni- 
formly track the if-band flux. 

For comparison, we also estimate the fraction of A"-band 
light enclosed by isophotes passing through the locations 
of 135 SNe la in early-type hosts discovered between Jan- 
uary 1990 and September 2011. The SN types, coordinates 
and host galaxy identifications were obtained from the IAU 
Central Bureau for Astronomical Telegrams (CBAT) web- 
sitqj. A SN la is included in the sample if its host galaxy 
has a morphology type of E, SO or SAO as classified in the 
NASA/IPAC Extragalactic Database (NED jf| and has a red- 
shift below 0.07. Our full SN la sample shows a deficit in the 
inner region of galaxies. To determine whether this reflects 
a true deficiency in the underlying distribution, we need to 
know the search efficiency as a function of projected galac- 
tocentric distance, especially for small distances. This factor 
is hard to determine precisely, since our sample is not drawn 
from a single survey with a well-understood selection func- 



2 http://www.cbat.eps.harvard.edu/lists/Supcrnovae.html 

3 http://ned.ipac.caltech.edu/ 
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Figure 1. Distribution of "calcium-rich" transients vs. enclosed 
fractional galaxy X-band light. The upper and lower limits are 
derived by fixing the galaxy Sersic index to 0.2 and 4 respectively 
(see text). The grey line shows the expectation if the transients 
follow the light distribution. The grey and dark dashed histograms 
represent respectively SNe la and SNe la found outside 20% K- 
band light. 



tion. We can proceed by making the reasonable assumption 
that the discovery efficiency drops towards the bright core of 
galaxy due to low contrast or blending, but is relatively flat 
outside 20% of stellar light. If we only include SNe la found 
outside 20% of stellar light (and assume 20% of events are 
cut by this criteria), the remaining SNe la track the light 
well, as seen in Figure [1] This agr e es wi th the result found 
in r-band bv iForster fc Schawinskil (|2008l), although Tf-ban d 
is a better tracer of mass than r-band l|Rix fc Riekdll993h . 
The decreased efficiency near galaxy centres has a minimal 
impact on our results regarding the "calcium-rich" sample. 
Using a di fferent fractional flux measurement, 
iKellv et all (120081 ) have shown that SNe la and SNe 
II follow the distribution of (/-band light while SN Ic and 
long G RBs are concentrated in the brightest regions of their 
hosts i|Fruchter et al.l I2006T) . The "calcium-rich" transients 
appear to be distributed differently with respect to host 
galaxy light from any of these samples. 



4 METALLICITY AND AGE CONSTRAINTS 

We now translate the distribution of the "calcium- rich" tran- 
sients with respect to host galaxy surface brightness into 
constraints on the age and metallicity of their progenitors. 
Galaxies are complex systems, often containing multi- 
ple stellar populations each with a different age and metal- 
licity. Central or average metallicities of host galaxies can 
be quite different from local metallicities at the supernova 
sites. Attempting to extrapolate the metallicity from the 
bright centre of a galaxy to its outskirts, where most of our 
events are located, is also dangerous since there is a larg e 
scatter in the gradient-mass relation l|Spolaor et al.l 120101 ). 
If the progenitors are very old, they may also have traveled 
far from the original formation site by the time of explosion, 
so the local metallicity may not give an accurate estimate 
of the progenitor metallicity. Finally, there is a significant 
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Table 1. Locations of calcium-rich transients in their host galaxies 



Transient 


Host 


Host 




Host 


Host r et ^ 


Distance from 


Fraction of Mass 


Fraction of Mass 


Name 




Typc^ 




M K 


(kpc) 


Host (kpc) 


Encloscd c 


Enclosed^ 


SN 2000da 


NGC 2768 


E 


-24. 8( 


-0.1, +0.2) 


6.2 


3.3 


0.58(-0.05, +0.11) 


0.43(-0.04, +0.08) 


SN 2001co 


NGC 5559 


SBb 


-24. 2( 


-0.1, +0.0) 


5.6 


6.6 


0.70(-0.08, +0.02) 


0.77(-0.08, +0.03) 


SN 2003H 


IC 2163 


Sc 


-24. 6( 


-0.1, +0.4) 


6.6 


6.7 


0.43(-0.04, +0.17) 


0.61(-0.06, +0.23) 


SN 2003dg 


UGC 6934 


So 


-23. 0( 


-0.2, +0.0) 


4.4 


1.4 


0.11(-0.02, +0.00) 


0.26(-0.05, +0.00) 


SN 2003dr 


NGC 5714 


Sc 


-22. 7( 


-0.1, +0.0) 


3.8 


2.5 


1.00(-0.10, +0.00) 


0.55(-0.06, +0.02) 


SN 2005E 


NGC 1032 


SO/a 


-24. 3( 


-0.2, +0.0) 


3.6 


23.9 


1.00(-0.04, +0.00) 


- 


SN 2007ke 


MCG +07-07-003 


E/S0 


-23. 1( 


-0.1, +0.0) 


0.9 


8.2 


1.00(-0.05, +0.00) 


- 


PTF09dav 


2MASX J22465295+2138221 


S? 


-23. 0( 


-0.3, +0.0) 


5.5 


40.7 


1.00(-0.03, +0.00) 


- 


PTFlOiuv 


CGCG 170-011 


s? 


-23. 4( 


-0.2, +0.0) 


2.6 


35.6 


1.00(-0.02, +0.00) 


- 


PTFllbij 


IC 3956 


s? 


-24. 1( 


-0.1, +0.0) 


3.1 


34.0 


1.00(-0.04, +0.00) 


- 



a Major axis of the ellipse enclosing half of the galaxy i^-band light. 
" Classification from HyperLeda llPaturel et al J - 20031. 
c Estimated in elliptical apertures defined by the transient locations. 
Estimated in circular apertures defined by the transient locations. 



scatter in the metallicity distribution function at a given ra- 
dius, i.e., there are some metal-poor stars in the centre and 
metal-rich stars in the outskirts of galaxies. To take all of 
these effects into account, we statistically compare the dis- 
tribution of projected galactocentric distances in our sample 
with the predicted distribution of stars with known metal- 
licity in self-consistent cosmological simulations. 

The simulations used in this comparison include the rel- 
evant physics such as hydrodynamics, star formation, super- 
nova feedback, a nd chemical enrichment , the details of which 
are described in lKobavashi et al.l (|2007l 1. Here we use a new 
simulation run with a resolution of 5 x 10 6 Mq, WMAP-5 
cosm ology and the Kroupa initial mass function (jKroupal 
l200ll l. A star formation parameter c = 0.02 and number 
of feedback neighbours (Nfb = 576) are used in order to 
match the observed cosmic star formation rate history and 
the mass-metallicity relations of galaxies. Metallicities, ele- 
mental abundances, ages, and positions of stars are recorded 
for analysis. Collections of particles in the simulation are 
identified as galaxies with the friend-of-friends method. 

We define metallicity as the logarithm of the mass frac- 
tion of elements heavier than helium compared to that of 
the Sun ([M/H]). From the simulations, we calculate the 
metallicity distributions as a function of projected enclosed 
stellar mass within host galaxies. We use circular apertures, 
since the mass distributions are approximately spherical in 
the outer regions of galaxies. We integrate the stellar masses 
along the line-of-sight since the observational data are mea- 
sured in projected radii. To be consistent with the simula- 
tions, we remeasure the enclosed Tf-band flux of our peculiar 
events in circular apertures. In cases where the transient is 
too far from the host to accurately determine the enclosed 
flux, we use the flux ratios from the previous section (see 
Table [TJ. We expect the choice of aperture shape to have a 
minimal effect on our results. 

The observed K-band total luminosities (Lk) of the 
host galaxies are converted to stellar masses for compari- 
son with simulations. Mj Lk values of 0.8, 0.7 and 0.5 are 
used for E/SO, SOa/b and Sbc/d type galaxies respectively 
l|Bell fc de Jong||200ll l and M/L K of 0.5 is used for the three 
hosts of unknown type. The resulting mass range is about 
1.2 x 10 10 M Q to 1.4 x 10 11 M , and therefore we select the 



galaxies in the simulation whose total stellar masses fall in 
the range 1 x 10 10 M Q to 1 x 10 11 M . Note that the result 
does not change qualitatively if a wider range of galaxies is 
included, as the stellar mass distribution peaks within the 
selected host mass range. 

Figure [5] shows the distribution of "calcium-rich" tran- 
sients compared to the stellar distribution for different 
metallicity ranges inferred from our simulations. In general, 
the "calcium-rich" transients lie further away from the host 
than stars of [M/H] <C —0.5, consistent with a relatively low 
metallicity origin. 

If only the five events in the K12 sample are considered, 
the distribution is skewed further from the host and is com- 
pletely different than any of our stellar metallicity curves in 
Figure [21 One way to resolve the discrepancy is to introduce 
an age limit on the stellar population. In the CDM picture, 
old stars form in small galaxies, which successively merge to 
form a large galaxy. Old stars that formed in a deep potential 
well tend to be located in the galactic bulge at present, and 
the chemical enrichment timescale is so short that the metal- 
licities of the majority of such old stars are not very low. On 
the other hand, old stars that formed in a shallow potential 
well, whose parent galaxy/cluster is disrupted, tend to be 
located in the outer halo at present, where no young stars 
form. 

Our sample does not have the statistical power to con- 
strain the age. For illustrative purpose, we set a cutoff age 
of 10 Gyr, corresponding to the lifetime of a star with zero 
main-sequence mass of about 1 M©. If we only trace stars 
that are older than 10 Gyr in our simulation, a dominant 
fraction of the metal poor population is found at large dis- 
tances from the galaxy centres. As can be seen in the bot- 
tom panel of Figure [31 even at this old age, the curve that 
comes close to the K12 sample has extremely low metallicity 
([M/H] £-3). 



5 CORRELATION WITH GLOBULAR 
CLUSTERS 

Another possibility leading to such large offset for classes 
of objects is if they are associated with globular clusters 
(GCs). GCs are detected throughout the extended halos 
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Figure 2. Observed distribution of the peculiar transients com- 
pared to those of stars of different metallicity within galaxies as 
predicted by our simulation. The full sample is shown in the top 
panel. In the bottom panel, the shaded area represents the con- 
tour for the K12 sample. Note that in the bottom panel, the 
x-axis range is zoomed in by a factor of 5 to enable the reader 
to distinguish visually between curves corresponding to different 
models. 



Figure 3. Observed distribution of the peculiar transients com- 
pared to those of stars of different metallicity and older than 10 
Gyr within galaxies as predicted by our simulation. The full sam- 
ple is shown in the top panel. In the bottom panel, the shaded 
area represents the contour for the K12 sample. As in Figure [2] 
in the bottom panel the x-axis range is zoomed in by a factor 
of 5 to enable the reader to distinguish visually between curves 
corresponding to different models. 



of galaxies (som etimes over 100 kpc away from the centre, 
iBrodie fc Strad er 2006). The high stellar densities in GCs 
lead to an enhanced rate of close binaries and unusual ob- 
jects, which could provide unique channels for exotic tran- 
sients. In this section, we test whether the progenitor pop- 
ulation could be correlated with GC systems. 

Direct photometric observations have failed to detect 
GCs of typical brightness at the sites of some of our 
nearby events. For SN 2000ds, Hubble Space Telescope 
imaging revealed no progenitor brighter than My ~ —5.5 
(|Van Dvk et al.ll2003l ). ruling out an association with a typ- 
ical GC at M v ~ -7. For SN 2005cz, a likely "calcium-rich" 
member that is not included in our statistical study, a non- 
detec tion down to Mv ~ —6.9 was obtained IjPerets et al.l 
I2011T ). This turns into a fairly restrictive limit, particu- 
larly considering that brighter GCs ten d to be more mas - 
sive and more centrally concentrated (cf. lMcLaughlinll2000r ) 
and thus more likely to harbour unusual objects. Because 
the calcium-rich class is still poorly defined, and there ex- 
ists potentially other large-offset transient populations (e.g., 



PTFlOops) we, despite these photometric limits, carry out 
the following comparison. This analysis should be useful for 
further studies where observations of comparable absolute 
depth are not achievable for more distant supernovae. The 
method described can also be generalised and applied to 
statistical studies of other transient classes which currently 
have too few members for study. 

For each host, we attempt to model the GC profile and 
calculate the fraction of GCs within (or outside) the tran- 
sient location, as we did for the stellar light in section [3] 
The projected radial distribution of globular clusters in a 
galaxy is often described by a power law, with s lopes be- 
tween ~ -2.5 and ~ -1.5 (|Brodie fc Stradeij[2006l . and ref- 
erences therein). More luminous galaxies tend to have shal- 
lower profi_les_while no obvious trend is found with galaxy 
type IJHarrislI 19861 ). This luminosity-slope relation is at least 
partly related to the fact that brighter galaxies are often 
larger in their physical extent. 

Since the GC radial distribution tends to flatten at the 
centre of a galaxy, we describe the GC projected surface 
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Figure 4. Cumulative distribution of "calcium-rich" transients 
as a function of enclosed fraction of GCs. The dashed line shows 
the expectation if the transients follow the GC distribution. The 
dotted line is for a fixed GC profile core size of 2 kpc and the 
dash-dotted line is for a cutoff of GC extent at 100 kpc. 



density at radius r as po(l + ( r / r c)~~ a )~ 1 ■ The core size 
r c is derived using the relation illustrated in Figure. 10 of 
iForbes et al. (119961 ). The slope a is estimated with Eq. (11) 
in lHarrisI ( 19861 ). Integrated V-band magnitudes of the 
galaxies are obtained from the NED database or converted 
from SDSS g and r model magnitudes. For three hosts with- 
out V^-band or SDSS photometry, we assume (V — K) — 3.2 
mag (the median for the other seven galaxies). We also ap- 
ply a realistic but somewhat arbitrary cutoff for the GC 
population extent of 120 kpc. The projected distance be- 
tween a transient and the centre of its apparent host is de- 
rived using the angular separation and the CMB rest-frame 
distance in the NED database, assuming Ho=70 km s _1 
Mpc -1 , fi m =0.27 and $!a=0.73. The resulting distribution 
is shown in Figure [4] While the histogram seems to suggest 
the "calcium-rich" transients may have a slightly flatter dis- 
tribution than the GCs at the inner regions of the galaxies, 
the discrepancy is not statistically significant. We show that 
the result is not very sensitive to the choice of cutoff dis- 
tance as all of our transients lie within 50 kpc. On the other 
hand, we see a stronger dependence on the core sizes. Given 
the many uncertainties in our GC profile estimate, we can 
only conclude that we do not rule out the possibility that 
"calcium-rich" transients originate in GCs. 

The number of GCs per unit stella r mass increases wit h 
decreasing galaxy mass at low masses (JMiller fc Lotal2007h . 
This tren d, combined with th e galaxy stellar mass function 
(e.g., from Bald rv et al.ll2012r ). means that a significant frac- 
tion of GCs exist in relatively faint/small galaxies. Unfor- 
tunately, seven of our 10 SNe were discovered in the LOSS 
survey, which targets large nearby galaxies a nd which may 
under-represent low surface brightness hosts (JLeaman et al.l 
I2011T ). The tight distribution of K-b&nd magnitudes of the 
identified hosts roughly agree with the typical luminosity of 
a LOSS galaxy. For a GC-related sample collected entirely 
through untargeted surveys in future, we would expect to 
see a wider range of host masses/luminosities. 



6 DISCUSSION 

6.1 Sample Bias 

Although the PTF "calcium-rich" transients all occurred far 
from the centres of their respective host galaxies, this is not 
a selection criterion for the class. The survey strategy can 
potentially introduce a bias, but PTF regularly discovers 
SNe at a range of galactocentric distances, suggesting that 
selection bias is not primarily responsible for the observed lo- 
cations of the "calcium-rich" transients. On the other hand, 
the targeted LOSS survey imposes a bias in their sample 
against objects lying far away from any bright host. By sim- 
ply combining events from PTF and LOSS, we run the risk 
of mixing two samples drawn from different populations. A 
larger sample drawn from ongoing and future untargeted 
surveys will improve the situation. For now, we try to verify 
that each of our sub-samples is statistically consistent with 
the overall population. 

The LOSS survey monitors a predefined nearby galaxy 
population with t he 0.76-m Katzman A utomatic Imaging 
Telescope (KAIT) j|Filippenko et al.ll200Jr ). We estimate the 
equivalent scale radii covered by KAIT's relatively small 
(6.7 x 6.7 arcmin 2 ) field of view at the distances of the 
seven LOSS discoveries. For the 2 (4) nearest events, this 
value is less than 30 (40) kpc. KAIT may therefore have 
missed transients with separations from their host galaxy 
larger than this. 

If the 3 events from PTF are drawn from the distribu- 
tion of all ten "calcium-rich" transients, we expect another 
1 (2) events detected inside about 50% (70%) stellar mass. 
Recognising a low luminosity "gap" transient near the bright 
core of a host galaxy can be challenging. For example in the 
host of PTFlOiuv (r = 14.2), a 19th magnitude object would 
be ~25% brighter than the underlying host at around 70% 
stellar mass. However, the number of transients potentially 
missed is still roughly consistent with fluctuations from Pois- 
son statistics given the small sample size. 



6.2 Implications 

A metallicity effect remains an appealing explanation for 
the large offsets of the "calcium-rich" transients. In either a 
thermonuclear or a core-collapse model, metallicity can af- 
fect the evolution of a progenitor star, impacting on its age, 
structure and chemical composition at the time of explosion. 
For example, production of a normal SN la (in a single- 
degenerate scenario) is inhibited at low metallicity, if the 
wind from the accreting white dwarf is not strong enough to 
stabilise the mass transfer (jKobavashi et al.lll998l ). Perhaps 
an exotic sub-luminous explosion can occur under these cir- 
cumstances. In the case of a massive star, wind-driven mass 
loss is closely related to metallicity. At lower metallicity, the 
mass loss becomes weaker and a larger helium core is re- 
tained, which may form a black hole dire ctly, resulting in 
a faint supernova (jKobavashi et al.l 1201 1! ) or no explosion 
JHeger et al.ll2003l) . 

Moreover, the extreme cases of "calcium-rich" tran- 
sients (the K12 sample) seem to lie further away from the 
host than the region where metallicity falls steeply. An ad- 
ditional parameter may be needed to interpret this distribu- 
tion. While we cannot strongly constrain the age with our 
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sample, our simulation show that extremely metal poor stars 
older than 10 Gyr are expected to be predominantly found 
at the same distant locations in galaxies as the K12 tran- 
sients. This trend also argues strongly against any massive 
star origin. This conclusion is consistent with the fact that 
some of the hosts are early-type galaxies without evidence 
for recent star formation. 

We also investigated if these peculiar transients could be 
associated with globular cluster systems, which often have a 
more extended spatial distribution than stellar light. In sec- 
tion[5]we showed that the radial distribution of the "calcium- 
rich" transients is consistent with the GC distributions. 
However, deep photometric observations argue otherwise 
and the discussed sample is not likely to originate from GCs. 
The GC connection may be more relevant for tran sients 
with suspecte d binary progenitors (|Pfahl et al.l 120091 ) . such 
as PTFlOops (|Maguire et alJuOllT) . One possible interpreta- 
tion for this unusual subluminous SN la is tha t it arises from 
the m erger of two equal-mass white dwarfs (|Pakmor et al.l 
[201CJ). Sufficiently violent mergers might preferentially oc- 
cur in the dynamical environment of globular clusters where 
non-sychronized and/or eccentric bina r ies co uld be more 
common. Alternatively, iRosswog et al.l (J2009J) have shown 
that the rate of collision-induced thermonuclear explosions 
is significantly enhanced in the dense cores of GCs. 

Finally, it is wo rth noting that the classification in 
iKasliwal et al.l (|2012l ) allows for a diversity in the photo- 
spheric phase. The spectra of PTF09dav show similarities 
to those of SNe la, but the rest of the sample are better 
matched to SNe Ib/c. In the nebular phase, differences ex- 
ist in the detection of oxygen emission. Nevertheless, their 
observed features imply comparable ejecta mass and energy 
release. If the "calcium-rich" transients are produced by one 
explosion channel, the model has to be able to account for 
this range of behaviour. 



7 CONCLUSIONS 

We have carried out a statistical study of the locations of 
"calcium-rich" transients within host galaxies to look for 
clues of progenitor properties. In particular, this class shows 
strong preference for large offsets from their apparent host 
galaxies, and do not follow the stellar mass distribution 
within galaxies as approximated by A'-band light. For these 
outlying transients, direct measurement of properties of the 
local stellar population is impossible, due to low surface 
brightness. Extrapolations to such large distances are also 
problematic, as significant uncertainties are involved. For 
our study, we choose to utilise self-consistent cosmological 
simulations. We show that as a group, the "calcium-rich" 
transient locations are consistent with relatively metal poor 
progenitors ([M/H] < -0.5). 

The sample selected in IKasliwal et al.l i|2012j ) includes 
the most extreme cases. It is not clear at the moment 
whether this is related to some selection bias, or sub-groups 
within the overall "calcium-rich" class. If these extreme 
cases are representative of a distinct class of explosions, their 
remote locations may imply a combination of low metallic- 
ity and old progenitor age. Events related to massive young 
stars are consequently disfavoured. 

The radial distribution of the "calcium-rich" transients 



is nominally consistent with the GC distributions. However, 
deep non-detection limits of the progenitors, in the few cases 
where they are available, disfavour an association of the pro- 
genitors with GCs. 

We expect ongoing and future surveys to generate a 
larger sample of transients for statistical studies like our 
own, including multiple instances of newly discovered rare 
types of transient events (e.g., PTFlOops). Our method may 
also be used to quantify the dependence of the properties of 
normal core-collapse SNe on the metallicity of their progen- 
itors. 
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